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Abstract: We discuss the distinct approaches for high density QCD 
(hdQCD) in the asymptotic regime of large values of parton density. We de- 
rive the AGL equation for running coupling constant and obtain the asymp- 
totic solution, demonstrating that the property of partial saturation of the 
solution of the AGL equation is not modified by the running of the cou- 
pling constant. We show that in this kinematical regime, the solution of the 
^ ■ AGL equation coincides with the solution of an evolution equation, obtained 

O recently using the McLerran-Venugopalan approach. Using the asymptotic 

O 

o 



> 

oo 



behavior of the gluon distribution we calculate the F 2 structure function 
assuming first that the leading twist relation between these two quantities 
is valid and second that this relation is modified by the higher twist terms 
associated to the unitarity corrections. In the first case we obtain that the 
corresponding F 2 structure function is linearly proportional to Ins, which 
agrees with the results obtained recently by Kovchegov using a distinct ap- 
proach. In the second softer behavior is obtained. In both cases, the 
^ ■ F 2 structure function unitarizes and the Froissart boundary is not violated 
in the asymptotic regime of high density QCD. We conclude that the partial 
saturation of the gluon distribution and the unitarization of the structure 
function are general results, well-established by three distinct approaches for 
high dense systems. 
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1 Introduction 



The parton high density regime in deep inelastic scattering is one of the 
frontiers in perturbative QCD (pQCD). This corresponds to the small x 
region and represents the challenge of studying the interface between the 
perturbative and nonperturbative QCD, with the peculiar feature that this 
transition is taken in a kinematical region where the strong coupling constant 
a s is small. By the domain of perturbative QCD we mean the region where 
the parton picture has been developed and the well-established methods of 
operator expansion and renormalization group equations have been applied 
sucessfully. The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions IJ, which are based upon the sum of QCD ladder diagrams, are the 
evolution equations in this kinematical region. In the limit of small values 
of x (< 10~ 2 ), on the other hand, one expects to see new features inside 
the nucleon: the density of gluons and quarks becomes very high and a new 
dynamical effect associated with the unitarity corrections is expected to stop 
the further growth of the structure functions. Ultimately, the physics in the 
region of high parton densities will be described by nonperturbative meth- 
ods, which is still waiting for a satisfactory solution in QCD. However, the 
transition from the moderate x region towards the small x limit may be very 
well be accessible in perturbation theory, and, hence, allows us to test the 
ideas about the onset of nonperturbative dynamics. In the QCD at dense 
systems (hdQCD) we still have a small parameter (the coupling constant a s ), 
which implies that we can start to approach this regime using the developed 
methods of perturbative QCD. 

The expectation of the transition for the high density regime can be un- 
derstood considering the physical picture of the deep inelastic scattering. In 
the infinite momentum frame (IMF) the virtual photon with virtuality Q 2 
measures the number density of charged partons having longitudinal momen- 
tum fraction x and transverse spatial size Ax t < 1/Q- When Q 2 is large, 
a s (Q 2 ) is small, so that the struck quark can be treated perturbatively. Also, 
when Q 2 is large the struck quark is small, Axt ~ l/Q, so that one can pic- 
ture the struck quark as being isolated, far away from similar quarks, in the 
proton. Thus, so long as the parton distributions are not large, the partons 
in a proton are dilute. However, if the parton distributions get large enough, 
which happens when x is very small, partons in the proton must begin to 
overlap. If there is a sufficient amount of parton overlap then a given parton 
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will not act as a free quantum over its lifetime but will interact strongly 
with the other partons in the proton, even though a s may still be in the 
perturbative regime. 

The behavior of the cross sections in the high energy limit (s — > oo) and 
fixed momentum transfer is expected to be described by the BFKL equa- 
tion ||. The simplest process where this equation applies is the high energy 
scattering between two heavy quark- ant iquark states, i.e. the onium-onium 
scattering. This process was studied in the dipole picture 0, where the 
heavy quark-antiquark pair and the soft gluons in the limit of large number 
of colors are viewed as a collection of color dipoles. However, one of the main 
characteristics of the BFKL equation is that it predicts very high density of 
partons in the small x region. The power behavior of BFKL cross section 
violates the Froissart boundary, which implies that the BFKL Pomeron de- 
scribes only the pre-asymptotic behavior at not very large energies and in 
order to find the true high energy asymptotics in perturbative QCD we need 
to unitarize the BFKL Pomeron. The understanding and analytical descrip- 
tion of the unitarity corrections and consequently of the high density QCD 
(hdQCD) is currently one great challenge. 

About 17 years ago, Gribov, Levin and Ryskin (GLR) 0] have started 
the description of the high density systems in the double logarithmic approx- 
imation (DLA) of perturbative QCD. They argued that the physical process 
of recombination of partons become important in the parton cascade at a 
large value of the parton density, and that these unitarity corrections could 
be expressed in a new evolution equation - the GLR equation. This equation 
considers the leading nonladder contributions: the multiladder diagrams, de- 
noted as fan diagrams. The main characteristics of this equation are that 
it truncates the series in the expansion in powers of the density at the first 
nonlinear term and predicts the saturation of the gluon distribution at very 
small x []; it predicts a critical line, separating the perturbative regime from 
the saturation regime; and it is only valid in the border of this critical line. 
Therefore, the GLR equation predicts the limit of its validity, and as all trun- 
cation is valid in a limited range. In the last decade, the solution |5|, [7| 
and possible generalizations || § of the GLR equation have been studied in 
great detail. At the moment, there are basically three distinct approaches 

1 We should be careful with this statement, since it is possible that the region of validity 
of the GLR equation ends before saturation is reached M . 
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for the dynamics at high densities: 



the Ayala, Gay Ducati and Levin (AGL) approach |jlQ |, which consid- 



ers the multiple pomeron exchange in the DLA limit of perturbative 
QCD, considering as basic degree of freedom the usual partons (quarks 
and gluons). The starting point of this approach is the proof of the 



Glauber formula in QCD [|TTJ, which considers only the interactions of 



the fastest partons with the target. In |10[ the authors have proposed 
a generalized evolution equation (the AGL equation) which takes into 
account the interaction of all partons in the partons cascade with the 
target and have analysed the solutions of this equation considering the 
semiclasssical method. The main conclusion was that the unitarity cor- 
rections strongly modify the behavior of the gluon distribution in the 
small x region; 



the McLerran-Venugopalan et al. (MV- JKLW) approach [|T!| [13[ , which 
is based on the effective Lagrangian formalism for the low x DIS and 
the Wilson renormalization group. The basic degree of freedom is the 



gluonic field. In |13| the authors have derived a general evolution equa- 
tion for the gluon distribution in the limit of large parton densities and 
leading logarithmic approximation considering a very large nucleus. In 
the general case this evolution equation is a functional equation, which 
does not allow to obtain analytical solutions. Recently, these authors 
have considered the DLA limit on their result [[14] and have shown 
that the evolution equation reduces to an equation with a functional 
form similar, but not identical, to the AGL equation. The solution of 
this equation was discussed in Ref. [H| considering the semiclassical 
method; 



the Kovchegov (K) approach ||16|| , which is based in the resumation of 



the multiple pomeron exchange in the leading logarithmic 1/x approx- 
imation [LLA(l/x)] using the dipole picture, where the basic degree 
of freedom are quark- ant iquark dipoles. The author has obtained an 
evolution equation for the interaction cross section of the qq pair with a 
nucleus considering the multiple scattering of the dipoles, which unita- 
rizes the BFKL Pomeron. The solution of this evolution equation was 
analysed by the author in Ref. [17] and independently in Ref. |TI| . 
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In contrast with the MV-JKLW and K approaches, in the AGL approach a 



comprehensive phenomeno logical analysis of the ep [IE, 20. 21 and eA 1 22, p3| 



processes exists. The analysis of the behavior of the distinct observables 
for the HERA kinematical region using the Glauber-Mueller formula was 
presented in Ref. [|l|J- In this kinematical region the solutions from the 
AGL equation and the Glauber-Mueller formula approximately coincide JII] 



The results from these analysis agree with the recent HERA data and allows 
to make some predictions which will be tested in a near future. Our main 
conclusion was that the unitarity corrections cannot be disregarded in the 
HERA kinematical region. 

All these approaches reproduce the small-x limit of the DGLAP evolution 
equations in the DLA limit of pQCD and the GLR equation as a first order 
unitarity correction. However, important questions still need to be answered: 
Since at high densities higher orders cannot be disregarded, what is the 
correct framework to treat these unitarity corrections ? Is there a common 
limit between these approaches ? 

Although the complete demonstration of the equivalence between these 
three approaches is still an open question, a first step was given recently in 



the Ref. [24]. In that paper we demonstrated the equivalence between the K 
equation and the AGL equation in the DLA limit, which motivates the use 
of the AGL approach, which considers the medium effects associated with 
the high density present at small values of x and/or nuclear collisions at high 
energies, as the dynamics at dense systems. Our goal in this paper is to give 
a second step in the demonstration of the equivalence between the three ap- 
proaches. Since in the general case the comparison between these approaches 
is a very difficult task, we will consider the asymptotic regime of very small 
values of x (very high densities) and demonstrate that in this regime there is 
an agreement in the predictions of these approaches, establishing a general 
behavior for the F 2 structure function in this kinematical region, which does 
not violate the Froissart boundary. 

The paper is organized as follows: In the next section we briefly review 
the derivation of the AGL equation and the asymptotic solution for fixed a s 



in order to fix notation and to present known results [10]. Later we derive the 



AGL equation for running a s and consider the asymptotic solution of this 
equation. We demonstrate that the main feature of the solution of the AGL 
equation, the partial saturation of the gluon distribution, is not modified by 
the inclusion of the running of the coupling constant. In the last subsection 
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we assume the evolution equation proposed by Jamal Jalilian Marian et al. in 



Ref. [14| and consider its asymptotic limit of high densities, showing that in 
this limit the predictions for the gluon distribution of this approach coincides 
with the AGL result. We will assume then the universality of the x and Q 2 
dependence of the gluon distribution (the partial saturation) in the asymp- 
totic regime of high density QCD and calculate the asymptotic behavior of 
the F 2 structure function in Section 3. We will consider two possibilities to 
calculate the structure function using as input the gluon distribution. First, 
we assume that the leading twist relation between F 2 and xG is not modified 
by the unitarity corrections. In this case we demonstrate that our results co- 



incide with the result obtained independently by Kovchegov in Ref. |l7j , and 



that F 2 does not violate the Froissart boundary in the asymptotic regime. 
This is a general behavior predicted by all approaches using the leading twist 
relation. Second, we consider that the relation between F 2 and xG is modi- 
fied by the higher twist terms associated with the unitarity corrections. We 
demonstrate that the behavior of the structure function is modified by the 
presence of these terms, but the Froissart boundary is not violated. Finally, 
in Section 4 we present our conclusions. 

We have tried to make this paper self-contained, and therefore have in- 



cluded some of the material already included in the early works 0, [14 



2 The AGL Equation 

2.1 Derivation and Asymptotic Solution for fixed a s 

We start considering the interaction between a virtual colorless hard probe 
and the nucleus via a gluon pair component of the virtual probe. The cross 
section for this process is written as 

fjG* A _ f l dz f fll^f { Q^ ru ^ z) \^ 99+ A {z y t) f (1) 



10 J TT 

where G* is the virtual colorless hard probe with virtuality Q 2 , r t is the 
transverse separation of the pair, z is the fraction of energy carried by the 
gluon and tyf is the wave function of the transverse polarized gluon in the 
virtual probe. Furthermore, a" +A (z,r 2 ) is the cross section of the interaction 
of the gg pair with the nucleus. To estimate the unitarity corrections we have 
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to take into account the rescatterings of the gluon pair inside the nucleus. 
Using the Glauber-Mueller approach, which considers the interations of the 
fastest partons with the target, we get 

a o: = f l dz ffn /^|* ? - W ^. ri ,,)P2|l-e- M ' (I '4 ,S, ''»] , (2) 



7T J 7T 

where x' = xj (z r\ Q 2 ) (x is the Bjorken variable), b t is the impact parameter, 

-4 

S(b t ) = (A/nRj^e r a is the gaussian profile function and aff = 1S 
the cross section of the interaction of the gg pair with the nucleons inside 
the nucleus. It was shown [25] that 

off = ^(3a s (i)/4)vr 2 r t 2 xG^(x,4) , (3) 

where xGn(x, -4) is the nucleon gluon distribution. 

The relation a G * A (x, Q 2 ) = (A7i 2 a s /Q 2 )xGA(x,Q 2 ) is valid for a virtual 
probe G* with virtuality Q 2 . Consequently, using the expression of the 
squared wavefunction we obtain that the Glauber-Mueller formula for the 
interaction of the gg pair with the nucleus is written as 

xG A (x,Q) = -j x —J ^jj—2[l-e . ]. (4) 



Q 



The AGL equation can be obtained directly from the above equation 
differentiating this formula with respect to y — Inl/x and e = In Q 2 /Aq Cd . 
Therefore the AGL equation is given by 

d 2 xG A {x,Q 2 ) 2Q 2 rd 2 b t 



dyde n 



2 



[±*[1- e-i^fcWM] , ( 5 ) 

J 7T 



where the dependence of cr^f in the virtuality of the virtual probe results from 
the derivative. The nonperturbative effects coming from the large distances 
are absorbed in the boundary and initial conditions. This equation is valid 
in the double logarithmic approximation (DLA). 

In our calculations we use the Gaussian parameterization for the profile 
function, which implies that the integration over bt is straithforward. More- 
over, in what follows we will consider the nucleon case (A = 1). Then, we 
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get 



° 2 - rGU - Q2) 2Q2R2 -{C + ln[K G (x,Q 2 )] + E l [ KG (x,Q 2 )}} , (6) 
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dyde n 

where C is the Euler constant, E\ is the exponential function and the function 
k g is defined by 

k g (x,Q 2 ) = 2Q2R2 xG ^ x ^ 2 "> ' ( 7 ) 

and represents the probability of gluon-gluon interaction inside the parton 
cascade. Using the above definition for k g we can rewrite the expression (^) 
in a more convenient form (for fixed a s ) 

d 2 k G (y,e) dk G (y,e) a s N c 2 2 

— k-k 1 a = {C + ln [k g (x,Q ) \+Ei [k g (x,Q ) } 

oyoe oy n 

= F{k g ). (8) 

In what follows we will present the asymptotic solution of the above equation 
[ |T0|j . Before some comments are necessary: 

• This equation matches the DLA limit of the DGLAP evolution equation 
in the limit of low parton densities (k — > 0); 

• In first order of unitarity corretions [0(k 2 )] the AGL equation matches 
the GLR equation; 

• The AGL equation sums all contributions of the order of K n absorbing 
them in the gluon distribution. 

From the analysis of the AGL equation we get that this equation describes 
the region of large k (large parton densities) where the use of GLR equation 
is not a good approximation. 

We expect that the unitarity corrections start to be important in kine- 
matical region where k g > 1. Assuming k g = 1 and using the definition 
of k g [Eq. (0)] we can approximately estimate the behavior of the gluon 
distribution in the region where the parton densities become large 

k g = \ ^ xG (x,Q 2 ) = ^P^ . (9) 



Therefore in a first approximation the gluon distribution saturates at small 
values of x. Below we will consider the asymptotic regime and demonstrate 
that this behavior is modified in the region where ka > 1 (very large densi- 
ties). 

Analysing the structure of the Eq. (§) we see that it has a solution which 
depends only on y. In |H| was shown that this solution is the asymptotic 
solution of the AGL equation. In this case we have 

{V ' = F(K a G symp ) , (10) 



dy 



with the solution 



k a G symp (y=yo) F{k' g 



y-yo- (ii) 



The determination of the behavior of the solution (|IT]) in the general case is 
a difficult task. However, in the limit of high parton densities we have that 
F(K G ymp ) — > In K, G ymp . Consequently, 

k a G avmP {y) dn' G a s N c t 

= (y - Vo) ■ 12 
C mp (y=») lnK' G tt vy yu; v ' 

Therefore at large values of densities and y ^> y the asymptotic solution is 
given by 

t asymp / \ &s-N c OL s N c fio\ 

k G w = yiny ~ y ■ (1 3 ) 

7T 7T 

We have checked that this solution is a good approximation at very small 
values of x [(9(10~ 8 )], where we expect that the density of partons will be 
very large, generating a high dense system. Moreover, as the energy behavior 
of the k function result from the y = Inl/x factor, we will assume in this 
paper that the asymptotic solution is given by (|13D and use this expression 
in the analysis that follows. 

Substituting the definition of k g [Eq. ([?])] in the above solution, we have 
that in asymptotic regime the behavior of the gluon distribution is given by 
(for fixed a s ) 

xG(x,Q")^^^ln(-). (14) 

67T Z X 
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Therefore, the gluon distribution do not saturate at small values of x, but 
is linearly proportional to Ins (s ~ 1/x). However, this behavior is softer 



than predicted by the DGLAP equation (xG oc exp[ylnl/x]) and the BFKL 
equation (xG oc x~ x , A > 0). We have that the gluon distribution present 
a partial saturation in its behavior. Moreover, we have that the gluon dis- 
tribution is directly dependent of the free parameter R, introduced by the 
profile function. In general this parameter is identified with the proton ra- 
dius. However, R is associated with the spatial gluon distribution within the 
proton, which may be smaller than the proton radius (see discussion in Ref. 

H). 

An identical dependence in Q 2 and x for the behavior of the gluon distri- 



bution in the asymptotic regime was obtained in Ref. ||18|| , where the solution 
of the K equation in the DLA limit was considered. This result is expected 
due to the equivalence between the K and AGL approaches in the DLA limit 



demonstrated in the Ref. 24 



In Ref. PD| Mueller has argued that the factor Inl/x in the gluon dis- 
tribution is associated with the one loop level of the calculations and that 
beyond of this level the distribution has the same form of the Weizacker- 
Williams gluon distribution, which is independent of the energy. However, 
this statement contrast with the studies of the parton evolution in the non- 



linear region [|H], [Lq| ; where the saturation is not present. The understanding 
of the difference between these results is an important open question. In this 
paper we assume (JTJ]) as the behavior of the gluon distribution in the asymp- 
totic regime. 

In the next section we consider the consequence of (0) in the behavior 
of the structure function. Before we will consider the implication of the 
use of running a s in the asymptotic solution of the AGL equation and the 
comparison between the predictions of the AGL and MV-JKLW approaches 
in the asymptotic regime. 



2.2 The AGL equation for running a s 

The solution ([14]) was obtained from the AGL equation [Eq. (§)] in the DLA 
limit of perturbative QCD and fixed a s . Our goal in this subsection is to 
shown that although the AGL equation is modified by the running of a s 
the general behavior of the gluon distribution, i. e. the partial saturation, is 
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independent of the approximation used. 

Since the QCD coupling constant is given by a s (Q 2 ) = 4:7i/((3 e), where 
Po = 11 — 2/3nf (rif is the number of flavors) and e = InQ 2 / Aq CD , the 
relation between the gluon distribution and the function k g can be expressed 
by 

xG(x,Q 2 ) = ^^^eK G (x,Q 2 ). (15) 

iV c 7T 47T 

Using this result in Eq. (||) we obtain that the AGL equation with running 
a s is given by 

92kG ^ e) ■ a +1 \ d -^A = N ^ m {C + ln[. G (x,Q 2 )] 



dyde \e J dy n 

+ E x [k g {x,Q 2 )]} 

= H{k g ) . (16) 

The AGL equation at fixed a s is a direct consequence of the above equation 
in the limit of large e. Moreover, this equation also has a solution which 
depends only on y. In this case we have 

{y ' e} - € -H(k g ) (17) 



dy 1 + e 



with the solution 



(2/-2/o). (18) 



Following the same steps used in the previous subsection we get that the 
asymptotic behavior of the gluon distribution is given by 

xG(x,Q 2 ) = -^^^ln(±). (19) 
1 + e 3n 2 x 

As expected at large values of e (Q 2 ) this solution reduces to the solution for 
fixed a s . The main difference is the behavior at small values of e, where the 
prefactor in (19) is important. However, the partial saturation of the gluon 
distribution is not modified by the running of a s . It agrees with the result 
that the unitarity corrections are expected to be relevant before of the next 
to leading order corrections |26|. 
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2.3 MV-JKLW approach in the asymptotic regime 



Recently an approach to the evolution of dense partonic systems within the 
framework of the Wilson renormalization group and the McLerran- Venugopalan 
approach was developed in Ref. 



This approach results in a nonlinear 
evolution equation for the generating functional of the color charge density 
correlators, which is valid to leading order in a s at densities which paramet- 
rically do not exceed l/a s . The equation is fairly complicated since it only 
requires ordering in longitudinal momenta during evolution and puts no con- 



straint in the ordering of transverse momenta (we refer |T^, |14]] for details). 
However, some special limits were considered: at low densities this equa- 
tion reduces to the BFKL equation and as consequence to the DLA limit 
of the DGLAP evolution equation at large values of Q 2 . In Ref. [D^] the 
authors have considered the DLA limit of the functional evolution equation 
obtaining the following equation for the impact parameter dependent gluon 
distribution xG(x, Q 2 , b t ) 



d 2 xG(x,Q 2 ,b t ) N C (N C -1) „ 



dyde 



where 



K(x,Q 2 ,b t 



Or 



2a s 



( xj> l — \ Ei ( — 



xG(x, Q 2 , b t 



(20) 



(21) 



jc(N c - 1)Q 2 

The general solution and the bt dependence of this evolution equation was 



discussed in Ref. 15 







At large values of densities (large re) we can approximate the Eq. (pO"D by 



d 2 xG(x,Q 2 ,b t ) N C {N C -1) 



dyde 



Q 2 



(22) 



To simplify our analysis we will assume that the b t dependence of the gluon 
distribution can be factorized xG(x, Q 2 , b t ) = xG(x, Q 2 )S{b t ) and consider a 
central collision (b t = 0), with S(0) = 1/nR 2 . Consequently, the evolution 
equation for the gluon distribution at large densities is given by 



d 2 xG(x,Q 2 ) _ N C (N C - 1)tt 
dyde 2 



R 2 Q 2 



(23) 
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The solution of this equation is straithforward 

xG (x, Q 2 ) = Nc(Nc ~ 1)n R 2 Q 2 In . (24) 



Therefore the solution of the Eq. (|2~0D has the same Q 2 and x dependence of 
the AGL equation in the asymptotic regime. The difference in the prefactors 
is a function of the distinct normalizations and approximations used in the 
two approaches. The demonstration of the equivalence between the two 
approaches in all kinematical region is still an open question. 

The main conclusion of our analysis of these two approaches is the uni- 
versality of the x and Q 2 dependence of the gluon distribution (the partial 
saturation) in the asymptotic regime of the high density QCD. In the next 
section we will consider the consequences of the partial saturation in the 
behavior of the structure function, what is experimentally measured. 



3 The Asymptotic Behavior of the Structure 
Function 

The common feature of the BFKL and DGLAP equations is the steep increase 
of the cross sections as x decreases. This steep increase cannot persist down to 
arbitrary low values of x since it violates a fundamental principle of quantum 
theory, i.e. the unitarity. In the context of relativistic quantum field theory 
of the strong interactions, unitarity implies that the cross section cannot 
increase with increasing energy s above log 2 s: the Froissart's theorem 
In this section we will calculate the structure function in the asymptotic 
regime using as input the solution of the AGL equation [Eq. (Ill)] which 
considers the unitarity corrections. 



Lets us start from the space-time picture of the ep processes |p8| . The 
deep inelastic scattering ep — > e+X is characterized by a large electron energy 
loss v (in the target rest frame) and an invariant momentum transfer q 2 = 
— Q 2 between the incoming and outgoing electron such that x = Q 2 jlm^y 
is fixed. The general features of the time development can be established 
using only Lorentz invariance and the uncertainty principle. The incoming 
physical electron state can, at a given instant of time, be expanded in terms 
of its (bare) Fock states 

|e > p hys= ipe\e > +ipe~,\ei > +■■■ ■ (25) 
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The amplitudes ipi depend on the kinematic variables describing the states 



\i >, and have the time dependence exp(-iEit), where Ei = J2i ymf + Pi is 
the energy of the state. The 'lifetime' tj ~ 1/0^ — E e ) of a Fock state \i > 
is given by the time interval after which the relative phase exp[— i(Ei — E e )] 
is significantly different from unity. If n > Ra the Fock state forms long 
before the electron arrives at the nucleon, and it lives long after its passage. 
New Fock states are not formed inside the nucleon. Therefore, the scattering 
inside the nucleon is diagonal in the Fock basis. If the state \i > contains 
particles with mass rrij, energy fraction Xj and transverse momentum p t j, 
we have that the transverse velocities Vtj = ptj/xjE e are small at large E e . 
Hence the impact parameters (transverse coordinates) of all particles are 
preserved. 

In terms of Fock states we then view the ep scattering as follows: the 
electron emits a photon (|e >— > \ej >) with E 1 = v and p\^ ~ Q 2 , after 
the photon splits into a qq (\ej >^ \eqq >) and typically travels a distance 
l c 1/rriNX, referred as the 'coherence length', before interacting in the 
nucleon. For small x, the photon converts to a quark pair at a large distance 
before it interacts to the target; for example, at the ep HERA collider, where 
one can study structure functions at x ~ 10~ 5 , the coherence length is as 
large as 10 4 fm, much larger than the nucleon radii. Consequently, the space- 
time picture of the DIS in the target rest frame can be viewed as the decay of 
the virtual photon at high energy (small x) into a quark- ant iquark pair long 
before the interaction with the target. The qq pair subsequently interacts 
with the target. In the small x region, where x <C <^r> ^ ne 99 P a i r crosses 
the target with fixed transverse distance r% between the quarks. It allows 
to factorize the total cross section between the wave function of the photon 
and the interaction cross section of the quark- ant iquark pair with the target. 
The photon wave function is calculable and the interaction cross section is 
modelled. Therefore, the structure function is given by [^9 



E 2 (x,Q 2 ) = JdzJ ^ll\^( z , rt )\ 2 a^(z,r t ) , (26) 



where 



n)\ 2 = z\{\z 2 + (1 - zfY KMtf + < K (er t ) 2 } , (27) 



i 
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a em is the electromagnetic coupling constant, e 2 = z(l — z)Q 2 + m 2 , rrii is 
the quark mass, n/ is the number of active flavors, e 2 is the square of the 
parton charge (in units of e), K ,i are the modified Bessel functions and z 
is the fraction of the photon's light-cone momentum carried by one of the 
quarks of the pair. In the leading log(l/x) approximation we can neglect the 
change of z during the interaction and describe the cross section a qq ~(z,A/r 2 ) 
as a function of the variable x. Considering only light quarks (i = u, d, s) 
F 2 can be expressed by p5 



= ^E^/*£? «"<*.'.)■ (28) 

u,d,s 



We have introduced a cutoff in the superior limit of the integration, in order 
to eliminate the long distance (non-perturbative) contribution in our calcu- 
lations. 

Using the expression of o qq [^5|] [Eq. |3|], we get 

F 2 (x,Q 2 ) = ^J^ d -^xG(x,Q 2 ). (29) 

This equation is a leading twist relation which will eventually break down 
when we consider the higher twist terms in the evolution. Below we discuss 
a generalization of the relation between the structure function and the gluon 
distribution and the implications in the asymptotic behavior. However, be- 



fore we will assume that the relation fl29|) is valid and will determine the 
behavior of the structure function in the asymptotic regime. 

Using the solution of the AGL equation in the asymptotic regime [Eq. 
(HID] as input in the Eq. (p9|) we get 



F 2 (x,Q 2 )^^-R 2 Q 2 ln(-) . (30) 

IT 6 \X/ 

We see that the partial saturation of the gluon distribution implies that the 
F 2 structure function does not saturate at small values of x, but is linearly 
proportional to Ins. The main point is that this behavior does not violate 
the Froissart boundary. 

An identical result for the behavior of the structure function was obtained 



by Kovchegov in Ref . |17| , using as input the solution of the evolution equa- 



tion proposed in [ HI . This evolution equation (the K equation) considers the 
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multiple pomeron exchanges in the leading logarithmic 1/x approximation 



[LLA(l/:r)], which implies the unitarization of BFKL Pomeron. In Ref. [24 
we have shown that the K equation is equivalent to the AGL equation in the 
DLA limit of perturbative QCD. The equivalence between the asymptotic 
behavior of the structure function obtained from the AGL approach (DLA) 
and the K approach [LLA(l/a;)] shows that the asymptotic behavior (28) 
is a general characteristic which is independent of the approximations used 
in the calculations. Furthermore, as the predictions of the AGL and MV- 
JKLW approaches for the gluon distribution are identical in the asymptotical 
regime, we have that the asymptotic behavior of the Fi structure function is 
well-established by the three distinct approaches. 

A comment is in order here. If we assume that the behavior of the gluon 
distribution is given by @, which is valid in the region where the density 
starts to becomes large, we then get that the structure function saturates at 
small values of x (F 2 oc R 2 Q 2 ). However, as shown in the previous discus- 
sion, this result is valid in the limited kinematical region where the higher 
orders in the parton density are not important. For a longer theoretical and 
phenomenological discussion about the parton saturation see for example the 



Refs. 30, 31 



We now consider that the relation between the structure function and the 
gluon distribution is modified by the particular type of higher twist terms 
associated with the unitarity corrections (See also |32|] for a similar calcula- 



tion in the infinite momentum frame). First we obtain the expression for the 
structure function considering these corrections and after we estimate the 
behavior of F 2 in the asymptotic regime. 

We estimated the unitarity corrections considering the s-channel unitarity 
constraint in the interaction cross section of the quark-antiquark pair with 



the target pi. In this case we the structure function is given by 



F 2 (x, Q 2 ) = -L £ e J ff / d\{l - e -i^<M/r?)s<fc) } . (31) 

Z7T , J — tt 7TT+ J 

u,d,s q2 1 

The use of the Gaussian parametrization for the nucleon profile function 

b 2 

S{b t ) = ^-e - "^, where R is a free parameter, as before, simplifies the 
calculations. We obtain that the F 2 structure function can be written as 
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F 2 (x,Q 2 ) = ^J2^ ff d -^{C + ln(K q ( X y t )) + £iK(*,r t 2 ))} , (32) 

A7V U,d,S J q" 



where K q = 4/9kg = (2a s /3i? 2 ) ixr\ xGn(x, -7). This equation allows to 
estimate the unitarity corrections to the structure function in the DLA limit. 
Expanding the equation (|3~2"D for small K q , the first term (Born term) will 
correspond to the usual DGLAP equation in the small x region [Eq. (|29|)] . 
In the asymptotic regime (large n q ) we obtain 

F ^Q 2 ) ~ S 5>? lf° d ^H^y t )) • (33) 

u,d,s q2 i 

Using the asymptotic solution of the AGL equation we can determine K q at 
large values of densities, and so 

r2 2 f& d 2 r t , [4a, . ,1 



Z7T 2 j ii 7rr? 



(34) 



Therefore, considering the contribution of the higher twist terms in the re- 
lation between the structure function and the gluon distribution, we predict 
the following asymptotic behavior for the structure function 



F 2 (*,Q 2 )^^f In 



— ln{~) 
3 x 



(35) 



We see that the inclusion of higher twist term implies a softer dependence of 
F2 with the energy than obtained using the leading twist relation. However, 
in both cases the structure function does not violate the Froissart boundary 
in the asymptotic regime of high density QCD. The demonstration of the 
behavior (^) using the other approaches for hdQCD is still an open question. 

From the results for the structure function in the asymptotic regime we 
can see that this regime is characterized by the identity 

dF 2 (x,Q 2 ) 2 

~^Q^ = F2Mh (36) 

which is an important signature of the asymptotic regime of high density 
QCD. This regime should be reached for the case of an interaction with 
nuclei at smaller parton densities than in a nucleon, since ka = As x k>n, 
where kn is given by the expression ([/]). 
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4 Conclusion 



In this paper we have analysed the asymptotic regime of high density QCD. 
We have shown that the partial saturation of the gluon distribution is a well- 
established result from the AGL equation for running and fixed a s and from 
the MV-JKLW approach in this kinematical regime. Using the asymptotic 
behavior of the gluon distribution we calculate the F 2 structure function 
assuming first that the leading twist relation between these two quantities 
is valid and second that this relation is modified by the higher twist terms 
associated to the unitarity corrections. In the first case we have obtained 
that the corresponding F 2 structure function is linearly proportional to In s, 
which agrees with the results obtained recently by Kovchegov using a distinct 
approach. In the second case a softer behavior is obtained. In both cases, the 
F 2 structure function unitarizes and the Froissart boundary is not violated 
in the asymptotic regime of high density QCD. We conclude that the partial 
saturation of the gluon distribution and the unitarization of the structure 
function are general results, established in three distinct approaches for high 
dense QCD systems. However, the demonstration of the equivalence between 
the distinct approaches in all kinematical regions and of the generality of the 
behavior of F 2 when the higher twist are considered are still important open 
questions. 
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